Scattering of short laser pulses from trapped fermions 
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We investigate the scattering of intense short laser pulses 
off trapped cold fermionic atoms. We discuss the sensitivity 
of the scattered light to the quantum statistics of the atoms. 
The temperature dependence of the scattered light spectrum 
is calculated. Comparisons are made with a system of classi- 
cal atoms who obey Maxwell-Boltzmann statistics. We find 
the total scattering increases as the fermions become cooler 
but eventually tails off at very low temperatures (far below the 
Fermi temperature) . At these low temperatures the fermionic 
degeneracy plays an important role in the scattering as it in- 
hibits spontaneous emission into occupied energy levels below 
the Fermi surface. We demonstrate temperature dependent 
qualitative changes in the differential and total spectrum can 
be utilized to probe quantum degeneracy of trapped Fermi gas 
when the total number of atoms are sufficiently large (> lO''). 
At smaller number of atoms, incoherent scattering dominates 
and it displays weak temperature dependence. 
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I. INTRODUCTION 

Rapid advances of trapping and cooling of alkali atoms 
[Q have led to the recent dramatic achievement of Bose- 
Einstein condensation (BEC) l|-|5|. Focusing these 
highly successful trapping and cooling methods onto 
fermions instead of bosons offers additional rich oppor- 
tunities for studying degenerate fermionic atomic gases. 
The description of these magnetically trapped Fermi 
gases presents additional simplification as the Pauli ex- 
clusion principle forbids the unsuppressed low energy s- 
wave collisions among atoms in the same hyperfine state. 
Thus these dilute trapped atoms behave very close to 
the ideal Fermi gas so ubiquitous in physics textbooks 
for decades. Recent highlights from several leading ex- 
perimental groups have indicated that we are at the edge 
of being able to explore quantum degenerate Fermi gas 
[g|,0 inside laboratories. 

Once such a degenerate gas is achieved, how do we per- 
form diagnostic measurements on its properties ? A stan- 
dard method in atomic physics is to probe the gas with 
light scattering. This is the problem to be addressed in 
our paper. These spectroscopic light scattering methods 
have already been suggested for BEC in both the weak 
and strong field scattering regimes ^,^ . These early the- 
oretical investigations had limited impact on BEC exper- 



imental observations so far, partly because the resonant 
photon-atom interaction becomes a complicated many 
body problem when a condensate is involved. The recent 
dramatic demonstration of low group velocity of light 
propagation inside a condensate [0, the Bragg scatter- 
ing experiment |ll[| , and the surprising observation of 
super-radiant collective spontaneous emission from MIT 
]l2[ all call for more detailed applications of quantum 
field theory of photons interacting with atoms. In the last 
several years, light scattering off Fermi degenerate atoms 
have already been discussed by several groups; some of 
these investigations focused on the case of a distribu- 
tion of ideal atoms obeying Fcrmi-Dirac statistics [[l3[ ; 
while others considered more exotic state for a Cooper 
paired fermionic condensate [ |l4| . One notable feature for 
a Fermi degenerate gas is that the Pauli exclusion prin- 
ciple blocks (inhibits) scattering events for atoms into 
already occupied states 111] . This paper aims at detailed 
quantitative information of light scattering from trapped 
Fermi gases within parameter regimes of current/future 
experiments, specifically results explicitly demonstrating 
temperature and quantum statistical effects are sought. 
In a previous paper involving two of us (LY and ML) 
an optical method for detection of the properties of BEC 
was proposed |16| . This detection involves the limiting 
case of scattering short but intense laser pulses from a 
system of cooled bosonic atoms in a trap. In particular, 
the case of laser pulses with areas of 2tiK was investi- 
gated. (The pulse area to be defined later is proportional 
to the integral of the slowly varying envelope of the elec- 
tric field multiplied by the atomic dipole moment). It 
was shown that such pulses mainly cause cyclic Rabi os- 
cillations for atoms in their excited and ground states. 
Thus to zero-th order involves no photon scattering (ex- 
cept the stimulated interaction with the driving field). 
However, atoms can spontaneously emit photons while 
in their excited states, therefore, the previous zero-th 
order picture involving coherent evolution of all atoms 
can be interrupted by spontaneous emissions of individ- 
ual atoms. Thus scattered light from these emissions can 
be collected and their properties reflect to a certain de- 
gree the properties of the trapped atoms. It was shown 
that for the case of BEC, above the critical tempera- 
ture, Tc, the coherent scattering is very weak and is pre- 
dominately in the forward direction due to phase match- 
ing effects. However, below the critical temperature the 
number of scattered photons increases dramatically and 



the coherent scattering occurs within a soUd angle deter- 
mined by the size of the condensate. For sufficiently short 
27r/'ir-pulses the system is preserved even below Tc so that 
this method can be potentially used as a non-destructive 
probe of BEC. 

This paper presents a continuation of our investiga- 
tions to the case of a system of trapped fermionic atoms. 
The ideal gas of a Fermi-Dirac distribution is considered 
in this paper. The static properties of harmonic trapped 
ideal fermionic atoms as considered by several groups pre- 
viously will be used as inputs n% . The more subtle case 
involving a weak attractive interactions between atoms 
could potentially develop into a BCS type Cooper paired 
condensate Q , and whose pulsed light scattering prop- 
erties [Q will be explored in a future publication. In the 
present study involving ideal non-interacting fermionic 
atoms, we do not expect to see a dramatic change in 
the spectrum as one cools the gas since no phase tran- 
sition occurs in the trapped gases even at zero temper- 
ature. However, with a finite number of atoms (say 1 
million) and cooled far below the Fermi temperature ||] 
one does expect that the quantum statistics should play a 
role in the spectrum. Indeed, we demonstrate that when 
number of atoms are sufficiently large such that coherent 
scattering dominates, spectrum of scattered light exhibits 
qualitative changes in its temperature dependence, a sig- 
nature potentially useful to probe quantum degeneracy 
of the system. 

The organization of the paper is as follows; first we 
review the formulation as was presented previously |g] 
for trapped bosons. In the second quantized form, the 
only difference lies at the commutation relations between 
atomic operators. We then calculate the spectra of scat- 
tered light taking into account the fermionic nature of the 
atomic operators. In Sect. II, we present the formulation 
of our model system. The near resonant laser pulse is as- 
sumed to be intense so that its interaction with trapped 
atoms is taken to be the zeroth order coherent process 
for which exact analytical solutions for the dynamics of 
atomic operators are given. In Sect. Ill, the interaction 
of weak scattered light with trapped atoms is treated 
quantum mechanically with perturbation method. The 
zeroth order solution of atomic operators from Sect. II 
is used to evaluate the dynamics of scattered photons. 
In Section IV, our new results for the total and differ- 
ential spectrums of scattered photons are then displayed 
numerically. A discussion section follows where both the 
angular and spectral distributions are presented. We also 
find the total number of atoms scattered for typical ex- 
perimental parameter sets. Finally we conclude with a 
summary of the essential physics learned from this in- 
vestigation. In the appendix we show how we can re- 
write the expression for the spectrum of the coherently 
scattered light in a form that is suitable for numerical 
calculation. 

As will become obvious by the introduction of many 
newly derived analytic formulae, the trapped fermion 
case is numerically far more difficult to calculate than 



the previous bosonic one. At low temperatures the 
fermionic system consists of many stacked energy levels 
|17[ whereas the bosonic system nicely condenses to a few 
of the lowest energy levels. Also the numerical methods 
for infinite series summations used previously for bosons 
are only applicable for the high temperature regime of a 
Fermi gas. New methods are therefore developed to cope 
with the additional difficulties of the fermionic system. 

The first few sections of this paper follows closely to 
the BEC situation. Initially the treatment is formally 
the same for both bosons and fcrmions, the only differ- 
ence being the commutation relations. However, these 
differences become significant later especially with the 
calculation of the spectrum. 



II. THE MODEL 

We consider a system consisting of N fermionic atoms 
confined in a trap interacting with light. Let us write 
down its Hamiltonian in the Fock representation and in 
the second quantized form S : 

n m 

-f ^ ^ / dkg{k) [rjfi^^{k)giai^ern- e^^ + h.c] 
where the Franck-Condon factors are 



(2) 



with R denoting the position operator of the atom (its 
nucleus). We have used atomic units (unless otherwise 
stated), the rotating wave and the dipole approxima- 
tion. The atomic annihilation and creation operators 
for the n-th state of the center of mass motion of the 
atoms in the trap are denoted by gfi and gl respectively. 
Since these operators are associated with the atoms in 
the ground electronic state, for the case of a spherically 
symmetric harmonic trap potential, n has three compo- 
nents {nx,ny,nz) and energy cui = ujtin^ + Uy + n^) 
where ujt is the trap frequency. The size of the trap is 
related to the size of the ground state of the trap po- 
tential, a = yJl/2MijJt i^ = 1). The atomic annihila- 
tion and creation operators in the excited state e^j and 
e^ raay experience different trap potential from that of 
the ground state. However, from our analysis of short 
pulse scattering, we find that the particular shape of the 
excited state potential is unimportant since atoms only 
spend a very short period of time in the excited state 
Q. The electronic transition occurs at the frequency wq- 
Since we are treating an s-state to a p-state transition the 
excited state operators are vectors e^ and el^. Annihila- 
tion and creation operators for photons of momentum k 



and linear polarization e^: (fi — 1,2) are denoted by aj: 

and at . All atomic operators obey standard fermionic 

anti-commutation relations. g{k) is a slowly varying cou- 
pling which is dependent on k. Its relation to the natural 

linewidth is 7 = (87r^fcQ/3c)u?(fco) , with fco — ^o/c. For 
notational convenience, we will suppress the indices g, e 
for the internal states. The convention being that the 



indices 



represents the center of mass states in the 



electronic ground potential whereas m, m' denotes the 
center of mass states in the excited state potentials. 

Note that the strong resonant atomic dipole-dipole in- 
teraction resulting from the exchange of transverse pho- 
tons is included in Eq. (|^). 

If the system is driven by a coherent laser pulse then 
we may neglect spontaneous emission effects during the 
pulse if the pulse is sufficiently short and intense. We can 
then safely substitute the electric field operator entering 
the interaction Hamiltonian in Eq. (|l|) by a c-number. 
The pulses we intend to use should have duration tl < 
300 ps or shorter, i.e. width jl = 1/tl ~ 3 x 10^ - 10" 
Hz. A first estimate shows that 7^ ^ 7 — 2.5 MHz, so 
that spontaneous emission may be neglected during the 
interaction time of the pulse with the atoms. The current 
estimate for single excitation spontaneous emission is a 
far better one than in the bosonic case as there is no Bose 
enhancement. However, the fermionic nature may come 
into play at very low temperatures where the effective 
spontaneous emission rate will be greatly reduced due to 
suppression by the Fermi sea of ground levels. Therefore 
it is a good approximation to assume that the effects of 
dissipative spontaneous emission and dispersive dipole- 
dipole interactions are small compared to the coherent 
driving laser during the interaction between the atoms 
and the laser pulse. We can then replace the product of 
the electric field operator {S'^-^^) and the absolute value 
of the electronic transition dipole moment (d) by 



d£(+) 



0" Z^ 



dk g(k, fj.)e 



ik-R—ickt 



(3) 



The envelope of the laser pulse is defined by g{k,ij,). H. 
is the peak Rabi frequency of the laser pulse. Using the 
assumption that the pulse is a plane wave packet moving 
in the fc/, direction with a central frequency lol and a 
linear polarization e^, we obtain 



2"' 



d£(+) ^ ^e-ir[7L(i - kL ■ i?M)]e^'"-«-''^-*. (4) 



The time dependent profile of the pulse, T{'-f]^t) is chosen 
to be real and we assume a gaussian shape with a peak 
equal to one at t = 0. 

In order for Eq. (0) to be valid, we need g{k, /i) to 
vary in momentum of the order of jl/c ~ 10 — 300 m^^. 
However, the Franck-Condon factors rifi^{k) vary on the 
scale of 5k, on the order of 1/a ~ 10^ m~^ for low n. 
For higher n's, 5k scales as l/\/n, so it becomes ~ 10'^ 



m ^ for the highest energy levels that are still available 
in the trap. Thus, we have 5k 3> 71,/c for all 7l, so we 
may validly replace k by k^ inside rjfi^ikL). Inserting 
this substitution into Eq. (||) and Eq. (Q), Eq. (||) can 
be written as 

n m 

+ -^^ilLt) [ exp(ic^Lt) Y, al^L ■ In + h.c] , (5) 

n 

where we have re- written the operators in terms of anni- 
hilation and creation operators of wave packets of excited 
states which originate from the n-th state of the ground 
state potential. 



Jn — / ^ Vnrhyf^L j^m- 



(6) 



These annihilation operators and their conjugate creation 
ones also obey the standard fermionic anti-commutation 

relations, i.e. {fi,fi,^} = 5ruv5qq', with q,q' = x,y,z 

enumerating the components of the vectors ffi and fi, . 
Their energies lo^ also vary very slowly for their cor- 
responding states and therefore for each of the wave 
packets /,-j, /- , their energy can be approximated by 
uji + ujo + k\/{2M). This assumes that the atomic 
wavepackets in the excited state potential will not ex- 
perience much coherent oscillation or diffusion (i.e. are 
in a sense frozen in shape) within the duration of the 
laser pulse (r^ <C 1/wt)- 

The Heisenberg equations that follow from the Hamil- 
tonian dsh now becomes linear. Thus at resonance, 
lol ~ loq + k\l{2M), and in the rotating frame in which 



9n{t) 
?L ■ h (t) 



-i{iAji+LiJL)t f^ 



fa, they become 



'■i^'^{lLt)eL ■ fn{t), 
-i^'^{lLt)9n{t). 



(7) 
(8) 



They can be easily solved analytically for any pulse en- 
velope 

9n{t) = 9n{~oo) cos [A{t)] 

-ieL-M~(x^) sin [A{t)], (9) 

cl ■ fait) = -ign{~oo) sin [A{t)] 

+eL-fRi~oo)cos[Ait)], (10) 

with the pulse area 



A{t)^^fjT{^Lt')dt' 



(11) 



Thus we see that each of the n levels of the ground state 
oscillator (when populated) creates an independent wave 
packet /a which is a superposition of the excited state 



wavefunctions. The population oscillates coherently be- 
tween the n-th ground state and the corresponding ex- 
cited state wave packet. We can view the behavior of the 
system as a set of independent two- level atoms coherently 
driven by the laser pulse. By using a pulse whose area is 
a multiple of 27r the system will be left in the same state 
after the duration of the pulse. Clearly, as \n\ increases, 
the approximations become less valid, but they should 
hold very well for the lowest 10^ available states of the 
ground state potential. 

The linear relations /s and e^j remain the same as be- 
fore for bosonic atoms |8| and their inverse can be easily 
derived using the sum rules 



n 



(12) 
(13) 



They will allow us to express solutions for gfi{t) and e,7i(i) 
as given above in Eqs. (pi) and dl^) (and their conjugates) 
uniquely in terms of S(fj(— oo), e^i(— oo), etc. 

Since spontaneous emission rate 7 from a dipole al- 
lowed excite states is non-zero, we do in practice have 
photons scattered from the atoms. The detailed discus- 
sions were given in our previous work |M|. Following the 
same perturbative treatment as developed there. One 
can work out the spectra of time dependent light scat- 
tering as given below. 



III. SPECTRUM OF THE SCATTERED LIGHT 

In order to calculate the spectrum of the scattered light 
we need to first determine the initial conditions for Eqs. 
(0) and (nG) at t = —00. We assume that initially the 
ground state energy levels were populated according to 
the Fermi-Dirac distribution (FDD) for non-interacting 
atoms in the harmonic well ||l^. We also neglect the 
much weaker p-wave collision interactions as they are en- 
ergetically suppressed in the temperature range of inter- 
ests here. 

The mean number of atoms in the n-th level at i = — cxd 
is therefore 



iV,7 



(.^ffff) 



-/3wii 



/(I 



-/3wfi 



), 



(14) 



where /3 = l/ksT, and z is the fugacity. The relation 
J2n ^n = ^ determines z as a function of /3, uJt and N. 
We can now calculate the spectrum of scattered pho- 
tons by using the Hamiltonian (|l|). We derive the Heisen- 
berg equation for the photon annihilation operator. 



«fcM " "*^^"fcM 



iQ{k)^glern ■ €i:^i]firn{k), (15) 



and its Hermitian conjugate for al . These equations 

are now solved perturbatively with respect to the atom- 
photon field coupling g{k). The formal solution of Eqs. 



n,rh 

xf di'e-*^'=(*-*%-'j^.e^(i')gt(t')- (16) 



The perturbative solution is then obtained similar to the 
case of bosons [B| but with care to take into account the 
anti-commutation relations between fermionic atomic op- 
erators. 

The total spectrum of scattered photons is defined as 



C{k,^,)^lim{al(t)a.{t)), 



(17) 



and can be divided into coherent and incoherent parts, 



C{k,fl) = Ccoh(fc,M) +Cin(fc,^). 



(18) 



The coherent part from (or: (t)), as in the single atom 
case, is proportional to the square modulus of the Fourier 
transform of the mean atomic polarization. The incoher- 
ent part, however, is due to the quantum fluctuations of 
the atomic polarization. Although in usual experiments, 
only the total spectrum can be measured, the division 
into coherent and incoherent parts is meaningful since 
they have significantly different angular characteristics, 
as we show latter. The total spectrum is 



C{k,fi)^\g{k)\'\eL 



^fc^i 



/ , ['7nimi(^)J Vn2m2\k) 

■ni,n2,m2 

Y^ rin[rnAkL)[Vn'^rn2{kL)]* 



ni ,mi ,712 ,m2 



n\ ,n' 



l''"2 

/OO /"OO 

dti / rfi2e~*('='=-"^)*ie^(='=~'^^)*^ 
- oc J —oc 

X (flih) ■ eL5ru(ii)ffi(t2)4(i2) • ei). (19) 

In the perturbative limit, we insert now the solutions of 
Eqs. (|) and (|l^ into Eq. (^. Here, the perturba- 
tion parameter is the coupling constant, g{k), of inter- 
actions leading to spontaneous emission. As the laser 
pulse is assumed intense, its corresponding interaction is 
much stronger, while the scattered light is weak so that it 
needs to be treated quantum mechanically. The pertur- 
bation approach consists of using zeroth order solutions 
for atomic operators in the Heisenberg picture to find the 
first order quantum correction to the photon operator. 
This method has also been applied in previous studies of 
scattering problems in bosonic systems B. 

At t = —00, the Heisenberg picture coincides with the 
Schrodinger picture, so that we omit in the following the 
explicit time dependence of the operators at t = —00. 
Since initially all atoms are in the ground electronic state, 
we obtain 

(/I ih) ■ eLgnAti)giit2)fn',it2) ■ el) 



X sin [A{ti)] cos [A{ti)] sin [^(^2)] cos [^(^2)] 

+ {gl^EL ■ fn.eL ■ fXgnO sin^ [A{h)] sin^ [^(ia)] , (20) 

where the expectation values (• • •) with respect to the 
FDD remains to be evaluated. Thus the only difference 
between the bosons and fermions so far is the evalua- 
tion of the expectation values according to the relevant 
statistics and initial conditions. The formal expressions 
are the same. 

The single atom spectrum can also be written as a sum 
of coherent and incoherent parts [pi , 



Ccohik,fj.) = 5'coh(^^) ViVa??mi(fc- /cl) 



(27) 



S{m) — 5coh(ro) + Sini'!^), 



(21) 



with VD = Suj/"fL — (ck — (-ol)/jl, and for a hyperbolic 
secant pulse l/cosh(7Lt) of area 2tt, 

Scohi^) = \p{ko)We^^ ■ CLfiTmy cosh' {Tim/ 2), (22) 
S'in(tn) = |p(/co)P(efc^, • eL)^7rn7^/sinh^(7rn7/2). 

Since the results are only weakly dependent on a particu- 
lar pulse shape we shall present results for the hyperbolic 
secant pulse only. 

Using the relationship 

^ [Vn'mik)] *mrn{kL) = ^n' (^L " k) , (23) 

m 

we obtain 

C{k,^i)^ ^ [qnin'^{k-kL)\*rin^n'^{k~kL) 

X [{gl'^gn,gl^gn'^)Scoh{T^) 



+ (gleL ■ fn,£L- /1^.9B^)5in(ro) 



Making use of the properties of the FDD we find 



(24) 



(25) 



and 



(gl'/L ■ fnieL ■ fl^gn'J = SH^n2Sn[n'^Nft'^. (26) 



The difference between bosons and fermions appears in 
the last term of Eq. (|2^) ; 1 — N for fermions and 1 + N 
for bosons Isl. In driving the above expression we used 
the fact that iV| = {{glignY) = Nn for fermions since 
the occupation of any given level is either or 1. This 
implies that the quantum dispersion of fermions is given 
by 5Nl = Nn{l - Nn). 

Inserting the above expressions in Eq. (p3), and per- 
forming tedious, but elementary calculations we finally 
obtain analytic expressions for the spectra. In particu- 
lar, the coherent part is 



identical in form as the bosonic system, but now with Nft 
representing the mean occupation number for a fermionic 
system. 

The incoherent part of the spectrum is 

Cin(fc,^) = 5coh(n7)^^AfH(l - Nfi,)\rififi,{k~ fc^)] 

n n' 

+ iV5i„(w). (28) 

Again the difference appears in the 1 — iV term on the first 
line of Eq. (|28| ) compared with l-|-iV in the case of bosons 
Q. We note that the incoherent spectrum ( |2^ ) consists 
of three parts coming from: i) quantum dispersion of the 

occupation numbers 5Nl = {{g^gnY) - (5,1 ff«)^ ii) P'^O" 
cesses of creation of the n-th wave packet accompanied 
by annihilation of the n'-th one for n ^ n\ and iii) the 
single atom incoherent spectrum. Obviously, both co- 
herent and incoherent spectra reflect quantum statistical 
properties of atoms since they depend on A^^-j's, which are 
described by the Fermi-Dirac distribution for our system 
of fermions. The first term in Eq. (Eq) for the incoherent 
spectrum depend explicitly on the statistical properties 
of atoms. 

The total number of emitted photons can be obtained 
by integrating the spectrum, 



^t,ot = 



Y,jdkC{k,ii), 



(29) 



which could also be divided into coherent and incoherent 
parts. By fixing the direction of k and integrating over 
the azimuthal angle Lp one can also define an angular 
distribution of photons dNtot{9) [and, correspondingly 
dN,,Y,{9) and dNi^{e)] 

= sm{e)de J2 [ dip f k'dk C(fc, fj.), (30) 



where 9 is the angle between k and fc^. We can also 
choose to define an integrated spectrum by fixing |fc|, 
and integrating the spectrum over the full solid angle. 



IV. NUMERICAL CALCULATION OF THE 
SPECTRUM 

The analytical expressions for both the coherent and 
incoherent spectra as obtained in Eqs. (^ and (^ are 
not directly applicable to straight numerical summations 
as they involve triple and six-fold sums respectively. In 
the previously work for trapped bosons, the analogous 
expressions were in terms of a power series of the fugac- 
ity z. A subsequent resummation of an auxiliary series 



provided a fast convergent numerical approach. How- 
ever, this same technique is only applicable when z < 1 
corresponding to a high temperature limit for the Fermi 
gas. We therefore needed a new method to calculate the 
spectrum at low temperatures when z ^ 1, which is ex- 
actly the region where the interesting effects of quantum 
statistics becomes important. 

The coherent spectrum as expressed earlier in Eq. (|27 
consists of a triple sum. Similar to the case for bosons 
|g|, when z < 1, we can write the spectrum as a power 
expansion of z 



Ccoh(fc,A^) — Sco\i{^) 






(1 - e^'^"')3 

2 



X exp [ - -(fc - kLf(? coth(//3wt/2)] 



(31) 



Thus we have transformed a triple sum into a single sum. 
Also for very high temperatures when z <C 1, the above 
single sum converges quickly. On the other hand, the 
quantum degenerate low temperature limit for fermions 
corresponds to z — > cx) when zero temperature is ap- 
proached. We thus use an alternative method where the 
original expression for the spectrum, Eq. (p7|), is written 
in terms of the Laguerre and the generalized Laguerre 
polynomials £„(•) and 'C™(.). The properties of these 
polynomials are then used to reduce the triple sum into 
a single sum over the generalized Laguerre polynomials. 
The new form is then (see Appendix A for details). 



Ccoh(fc,/i) = 5'coh(tn)e"(*^-*'-)'''' 

CSO 



n=0 



where 



P(n) 



ze 



-fSuJtn 



1 + ze-^'^*" 



(32) 



(33) 



is the mean occupation number of fermions in any of the 
degenerate energy states with principle quantum number 
n. 

By using this new expression we achieve enormous sav- 
ings in computational time over the initial triple sum. 
However the added complexity of generating the gener- 
alized Laguerre polynomials implies that this method is 
not as fast as the power expansion method which was 
only valid for z < 1 [g[ . 

The evaluation of the incoherent spectrum is a more 
difficult numerical problem as we see that the second 
term of Eq. ( |2q ) involves a six-sum. This is due to the 
fact that incoherent photon emissions corresponds to dif- 
ferent final and initial motional energy levels (unlike the 
coherent case where the two are the same). For the case 
of z < 1 we can again use the power expansion approach 
to write the incoherent spectrum as 



— <S'coh(w) 



(„^yi+;2 



^ [1 _ e-('i+'2)/3'^t]3 



h,l 



X exp[- {k-kL)^a^f{(3,h,l2)], 



with 



fiPJi^h) 



(l-e-'i^'^')(l 



-l20l^t\ 



1 _ g-{h+h)0'^t 



(34) 



(35) 



We note that the minus '-' sign in front of the third term 
in Eq. (|34| ) is due to quantum statistics. In the bosonic 
case considered earlier m a plus '-I-' sign was involved. 
This approach helped to reduce the six-sum into a double 
sum. Again this limit is only applicable for high temper- 
atures (when average energy per atom is much greater 
than the Fermi temperature). In the general case, some 
alternative simplifications are possible. We rewrite Eq. 
@as 

- Scoh{rAj)^^NfiNfi'\f]Hn'(k- fci)] 

n n' 

= 7V5coh(tn)-|-iVS'in(ro) 

- S'coh(t^) ^ ^ 



X P\nc{nx +ny-\-nz,mj,-^ lUy 

I 1 2 

V T T T 



m^) 



with 



Pinc{n,m 



2g-/3wtn ^^—pbJtm 



and 



(1 H- ze-^'^*")(l + ze-'^'^*™) ' 



(36) 



(37) 



(38) 



Where we denote the j-th component oi Sk ^ k — k^ by 
Skj . The number of sums can be reduced by exploiting 
the symmetry of the scattering geometry. Because the 
scattering is symmetric about the axis along the incoming 
laser direction k^. We can set one of the 6k components 
in the plane perpendicular to the laser axis to be zero, 
i.e. choosing the laser to be aligned with the z-axis we 
can set ky = 0. The incoherent spect rum with the aid of 



this symmetry and the identity Eq. (A3) is then 

Ci„(fc,M) = NS,ohi^) + NSini^) 
- 5eoh(n7)e~l^'=l'"' 
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with 



Vin, m) = ^ Pincin + y,m + y). 



(40) 



This finally reduces the six-sum to a four-sum. 



V. DISCUSSION OF THE SPECTRUM 

In this section we discuss three methods for analyzing 
the spectrum. 1) we look at the part of the spectrum that 
reveals quantum statistics. This represents all the collec- 
tive effects on the spectrum. 2) we integrate over either 
the frequency (angle) to obtain an angular (frequency) 
spectrum respectively. 3) the total number of scattered 
photons is calculated to illustrate the overall scattering 
behavior as a function of the temperature and hence the 
degeneracy of trapped fermions. 



A. Form functions 

A useful part of the spectrum is the component that 
describes the quantum statistics as opposed to the single 
atom component of the spectrum. We call this compo- 
nent the "form function" as it is indeed related to the 
Fourier transform of the average density profile of the 
degenerate gas. Temperature dependent behavior of the 
scattering spectrum, the transition from classical statis- 
tics at high temperatures to the Fermi-Dirac one at low 
temperatures, will manifest itself in the form function. 

In Eq. (pTl), the quantum statistical component is the 
expression to the right of the single atom term Scohi'^)- 
Thus for z < 1 we have the following form function 
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X exp [ (fc - kLfa^ coth(m/3cjt/2)] 



(41) 



For atoms obeying Maxwell-Boltzman distribution, the 
above summation can be analytically evaluated to give 



T: 



^o^{5uj,e) = iV^exp -{k-kLfa^coih{(iL0t/2) 



(42) 



Thus, at high temperatures coherent spectrum of scat- 
tered light from fermionic atoms should coincide with 
the classical results. For the general case of arbitrary 
temperature we can rewrite Eq. (p2[) as 
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Figures [l[g represent the coherent form functions for 
temperatures of ksT/Ep = 1.36, and 0.0016 with a sys- 
tem consisting of one million atoms. Results are pre- 
sented as rescaled by 7V^, since T^^^{0,0) = iV^, inde- 
pendent of temperature. In all numerical simulations we 
have set the dimensionless parameter k^a — 12.5. For a 
resonant transition wavelength of ~ 800 (nm) , our choice 
corresponds to magnetic trapping frequencies of about 
27r X 300 (Hz), 2n x 50 (Hz), and 2n x 23 (Hz) for ^Li, 
"'"K, and ®^Rb respectively. The ground state trap size 
is ~ 1.6 (/im). We see that the form function becomes 
broader as the temperature drops. At high temperatures 
phase matching effects become dominant where destruc- 
tive interference attenuates the scattering except for a 
narrow cone in the forward direction. As ksT drops be- 
low the Fermi energy Ep, phase matching becomes less 
important and the role of the quantum statistics of the 
gas becomes increasingly significant. As a rule of thumb, 
the effects of quantum statistics comes into play when 
ksT is less than one-half of iJ/p, consistent with the effects 
seen in evaporative cooling [|l9[ and recent experimental 
studies M. We note that once the quantum statistic ef- 
fects become significant, lowering the temperature does 
little to the shape of the form function as the lowest en- 
ergy levels become close to being "stacked" , and start to 
block further filling due to the Pauli exclusion principle. 
Any additional cooling will therefore only result in filling 
of the abundant higher energy levels to enlarge the Fermi 
sea. 

The form function for the incoherent spectrum when 
z < 1 is retrieved from Eq. (p4) 
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xexp[-(fc-fci)V/(/?,Zi,;2)]. (44) 

The corresponding form function for the general case 
from Eq. (^9|) would be 
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Figures Hjj display the incoherent form functions for two 
different temperatures kBT/Ep — 1.36 and 0.6 for one 
million atoms. The form function at the cooler temper- 
ature (Fig. 0) is narrower along the frequency axis and 
have a higher sharper peak than that of Fig. |[ At the 
high temperature limit we can estimate the height of the 
peak in Fig. by using Maxwell-Boltzmann statistics to 
find that J^f„(0,0)/7V = N/{2kTf. Putting in the num- 



(43) bers for Fig. |; A^ = 10^ and l/{kT) = 4.036 x lO^^ we 
obtain J^l,{0, 0)/N = 8.2 x 10"^ which is consistent with 



the numerical peak of 8.0 x 10^"^. We also know that 
in opposing limit of zero temperature the form function 
has a value of N at the peak. At frequency regions far 
away from the central peak location, the form function is 
mainly governed by the rapidly decaying exponential fac- 
tor exp[— (/c 



kL)'^a'^] in Eqs. (43) and 



i§. This is so 

because the weighting factors V of the equilibrium atom 
number occupations more than off-set the asymptotically 
growing (x"/n!) behaviors of the Laguerre polynomials 
L" (x) . Consequently, higher order Laguerre polynomials 
don't play any significant role in determining the asymp- 
totic character of the form function, and the form func- 
tion simply decays exponentially which physically corre- 
sponds to suppressed scatterings into too wide frequency 
separations. On the other hand, angular behavior of form 
function shows periodical character. 

We see that the form functions are asymmetric along 
the frequency domain but remains symmetric in angular 
one. This angular symmetry is basically due to our previ- 
ous assumption on the cylindrical symmetry of the scat- 
tering, since the exponential factor exp[— (A; — fcL)^a^] de- 
pends on angle 6 through the even function of cos 9. The 
asymmetry along the frequency direction is also consis- 
tent with the asymmetry of exponential factor with re- 
spect to frequency a; = c\k\. 



B. Angular and frequency spectrum 

The angular and frequency spectra are of interest as 
they can be easily observed in an experiment. We have 
already defined an angular spectrum of photons dNtot (^) 
previously by Eq. pO|). Note that the domains of the 
variables; the azimuthal angle ip, polar angle and k the 
magnitude of k are: ip G [0, 2tt], 6 G [0, tt] and k e [0, oo]. 
It is desirable to look at the spectrum of the coherent 
and incoherent components of the spectrum separately 
as their behavior are significantly different. We do not 
explicitly write down the expressions for the coherent sit- 
uation as it is trivial; it involves simply an integration and 
sum over the polarizations of the product of the single 
atom contribution 6*0011 and the coherent form function 
^^qJj. The incoherent situation is, on the other hand, a 
little more complicated. Using the expressions for Cin 
from Eq. (gj) and Eq. ( |39| ) we have the following general 
expression 



dNiniO) = 7r(l + cos^ 6) sinOde 

X {N / [5coh(n7) -t- 5i„(tn)] k^dk 
Jo 

) 

S,oh{r:o)J^L{0,k)k^dk}, (46) 



To obtain the analogous expression for the angle inte- 
grated frequency spectrum one simply integrate over 9 
instead of k in Eq. ( ^ ) . 



Figures ^ and |6| show coherent angular and frequency 
spectra for one million atoms with a laser pulse width 
of Tz, = 10 ps. As shown in Fig. |^, the coherent an- 
gular spectrum is narrow with a range from about —2 
to 2 degrees. The solid, dashed, and dash-dotted curves 
corresponds to temperatures of ksT/EF = 0.001, 0.5, 
and 1.0 respectively. Curves for lower temperatures are 
broader and greater in magnitude than the high temper- 
ature dash-dotted curve. Figure ^ shows the frequency 
spectrum for the same three temperatures using the same 
curve formats. The overall features for the coherent spec- 
tra are the increase in magnitude and broadening of the 
scattering as the temperature drops. As explained pre- 
viously in the form function section, the effects of phase 
matching becomes diminished and quantum statistical 
effects becomes more dominant as the temperature cools 
leading to an observable broadening and an overall in- 
creased scattering. 

Finally Figs. and || show incoherent angular and fre- 
quency spectra for one million atoms with a laser pulse 
width of Ti = 10 ps. The incoherent angular spec- 
trum is far broader than the coherent one with a range 
from about —180 to 180 degrees, covering the complete 
polar range. We used the same curve formats as be- 
fore for the three temperatures. We see that all three 
curves are almost on top of each other for angular ranges 
90 < 6* < 180 and -180 < 6* < -90, which corre- 
sponds to back-scattering regime. The single atom spec- 
tra S'coh, 5'in are appreciable only within a frequency scale 
~ 10^7, much shorter than the frequency scale ~ 10^7 
for any significant change in form functions F^^^i and 
-F;^. Thus, frequency dependence of form functions can- 
not play a role in determining the behavior of the spec- 
trum of the scattered photons. We approximate from 
Eq.(p6|) that angular behavior of differential spectrum 
follows (1 + cos^ 9) I sin 6*1 [N - Fj^ (0, 9)] . Therefore, we 
conclude that the decrease in the first peak is due to the 
increase of F;^ at large number of atoms and at low tem- 
peratures. The second peak is not affected since within 
the frequency scale of single atom spectra, the form func- 
tion is already diminishingly small at those angles. The 
minus sign in A^ — Fj^(0, 0) is due to anti-commutation 
of fermionic atomic operators, thus is purely quantum 
statistical. In Fig. we see that again all three curves 
follow one another closely. However they now possess a 
dip at the resonance frequency. Because these spectrum 
consists of incoherent processes, unlike the case of coher- 
ent scattering at high temperatures, there is no phase 
matching effects. The dip in the frequency spectrum is 
mathematically due to the requirement that all the three 
temperature curves meet at to — loq equal zero. For this 
particular point the form functions and hence the spec- 
tra does not depend on the temperature, reminiscent of 
some kind of optical theorem H . 



C. Number of scattered photons 

The total number of scattered photons is also directly 
observable in an experiment. By calculating the total 
number of scattered photons from a simple pulse excita- 
tion we can study its dependence on the temperature. We 
again separate these into coherent and incoherent com- 
ponents. The number of scattered photons scales as iV^ 
for the coherent case and as N for the incoherent one. 
For larger numbers of atoms one expects that the scat- 
tering will be dominated by the coherent scattering in 
the short pulse limit. Similarly, for low atom numbers, 
the incoherent scattering will become dominate and thus 
be observable when only few atoms are scattered. 

The coherent scattering becomes more dominant at 
lower temperatures for the system of one million atoms. 
Our preference for a dominant coherent scattering is due 
to the more sensitive nature of the coherent scattering 
to changes in temperature. This is shown clearly in the 
following figure, Fig. H showing coherent and incoherent 
scattering as a function of the temperature for one mil- 
lion atoms. The coherent results are plotted as circles 
while squares denote the incoherent ones. For compari- 
son the triangle plots are calculated using the Maxwell- 
Boltzmann distribution (MED) for a classical gas. The 
coherent curve displays temperature-dependent sensitiv- 
ity across the full range of temperatures both above and 
below Ep. The incoherent curve, on the other hand, is 
flat across the entire range of temperatures that we used. 
As is expected, the MBD and FDD coherent curves are 
the same at high temperatures since quantum statistical 
nature of the atom becomes less important as ksT/Ep is 
greater than 0.5. They start to deviate from each other 
for ksT/EF between and 0.5. This becomes very pro- 
nounced at very low temperatures when ksT/Ep is less 
than one-tenth. The FDD curve (circles) displays a flat- 
tening off at very low temperatures. The MBD curve, on 
the other hand, does not share this feature, but it dis- 
plays a rapid increase in scattering at these low temper- 
atures. Thus the flattening off is caused by the fermionic 
nature of the atoms. Previous work on a bosonic system 
displayed a dramatic increase scattering at these low tem- 
peratures. This flattening originates from the inhibition 
of spontaneous emission far below the Fermi tempera- 
tures where the lower levels are "stacked" so that spon- 
taneous emission into these levels are forbidden by the 
Pauli's exclusion principle fla]. 



theory, 2t:K pulses leave the system of trapped atoms 
practically unperturbed. At high temperatures when the 
average energy per atom is many times the Fermi temper- 
ature Tp the coherent scattering is very weak and occurs 
in a very narrow cone in the forward direction due to 
phase matching effects, which is the same effect as for 
the case of bosons m. As the temperature ksT becomes 
of order less than one half of Ep, angular distributions of 
scattered coherent photons broadens as the influence of 
the phase matching effects are reduced. Incoherent scat- 
tering is insensitive, in both angular and frequency spec- 
tra (hence total photons scattered) to the temperature 
at small atom numbers (< 10^). For such small number 
of atoms incoherent scattering dominates and scattering 
of short laser pulses is weakly temperature dependent. 
On the other hand, we find for larger number of atoms 
(> 10^ atoms with other parameters chosen as in this 
study), the angular incoherent spectrum, as well as the 
total number of scattered photons, exhibit qualitative ef- 
fects depending on temperature. 

The number of scattered photons increases as the sys- 
tem of atoms are cooled but the rate of this increase tails 
off as the fermionic nature sets in to suppress photon 
scattering events leading to atoms into already occupied 
motional states. This temperature dependent property 
is compared with results calculated for atoms obeying 
the Maxwell-Boltzmann statistics. In this case of distin- 
guishable classical atoms, the number of scattered pho- 
tons follow the fermionic system at high temperatures, 
as it should, but starts to deviate at low temperatures 
(around one half of Ep) where it continues to increase at 
a higher rate. Scattering of short laser pulses on a system 
of trapped atoms thus provides a useful method for de- 
tecting the temperature and hence degree of degeneracy 
of a fermionic system. 
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VI. CONCLUSIONS 

We have studied theoretically the scattering of intense 
short laser pulses off a system of cold atoms. We have 
presented a detailed theory of such processes. We have 
demonstrated that by scattering pulses of area 27tK one 
may observe signatures of fermionic degeneracy in a sys- 
tem of trapped atoms. In the regime of validity of our 



APPENDIX A: COHERENT SPECTRUM 

Starting from Eq. ( p7| ) we explicitly write down the 
triple sum expression for the coherent spectrum as 



Ccoh(fc, M) ~ Scoh{'^)\ 
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Since Rj is the position operator in the j'-component 
(j = x,y,z), the above matrix element is simply 
the diagonal elements of the displacement operator, 
{nj\'D{~iSkja)\nj) . We can then write this matrix el- 
ement in terms of the Laguerre polynomials using the 
following identity [pO[, 

{n\Vmm) = J^e-^i^'/'C-'C^-^Qe) (A3) 
V n\ 

to obtain 

{n,\Vi(,)\n,)=e-\^\^-/^Cn,m'), (A4) 



with ^j =^ —iSkja. We can now re-write Eq. (Al) as 
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where £, = {£,x,£,y,£,z)- We now rearrange the order of 
summation to obtain 



Ccoh(fc,M) = '5'coh(tn)e '^M ^ ^ 
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Since for the spherical symmetric trap under consider- 
ation, the mean occupation number Nn^^n ,n^ only de- 
pends on the sum of n^, riy and n^ (n — n^ + riy -\- riz), 
we can take it out of the inner sum and express it as 



P{n) 



'IJuJtn 



1 -I- ze-/^"'" ■ 



(A7) 



We then have 



Ccoh(^, ^J) = 5coh(tn)e-l«l' I Y, P[n) Y. 

X CnM.?)i^nMy?)^nA\S,z?)^ ■ (A8) 

Using the summation theorem of Laguerre polynomials 
PH to replace the inner sum with a single generalized 
Laguerre polynomial we finally obtain for the coherent 
spectrum. 



Ccoh(fc,Ai) - 5cohM| 5^F(n)£2(|^|2) 
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FIG. 1. The coherent form function for one million trapped 
fermionic atoms at a temperature of ksT /Ef = 1.36. 
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FIG. 2. The same as Fig. 1 but at a temperature of 
ksT/EF = 0.0016. 
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FIG. 3. The incoherent form function for one million 
trapped fermionic atoms at a temperature of ksT/EF ~ 1.36. 
Note that we have scaled the form function by A'' instead of 
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FIG. 5. The differential coherent scattering versus an- 
gle for ten thousand trapped fermionic atoms. The solid, 
dashed, and dash-dotted curves represent temperatures of 
ksT/EF = 0.001, 0.5 and 1.0 respectively. 
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FIG. 4. The same as Fig. 4 but at a temperature of 
keT/Ep = 0.6. 
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FIG. 6. The differential coherent scattering as a function of 
frequency for ten thousand trapped fermionic atoms. Follows 
the same curve format and temperatures as Fig. U. 
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FIG. 7. The same as Fig. ^ but for differential incoherent 
scattering. 
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FIG. 8. The same as Fig. M but for differential incoherent 
scattering. 
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FIG. 9. The total number of coherent and incoherent scat- 
tered photons as a function of temperature for one million 
trapped fermionic atoms. Coherent scattering calculated for 
a Maxwell-Boltzmann system is plotted as triangles. A 10 ps 
laser pulse width was used in the calculations. 



